arXiv: 1502.0495 lvl [nlin.AO] 17 Feb 2015 


cjji uiau 


Design and control of noise-induced synchronization patterns 

W. Kurebayashi 1 ^), T. Ishii 1 , M. Hasegawa 2 and H. Nakao 1 

1 Graduate School of Information Science and Engineering, Tokyo Institute of Technology - 2-12-1 O-okayama, 
Meguro-ku, Tokyo 152-8552, Japan 

2 Department of Electrical Engineering, Tokyo University of Science - 6-3-1 Niijuku, Katsushika-ku, Tokyo 125-8585, 
Japan 


PACS nn.mm.xx - First pacs description 
PACS nn.mm.xx - Second pacs description 
PACS nn.mm.xx - Third pacs description 

Abstract - We propose a method for controlling synchronization patterns of limit-cycle oscillators 
by common noisy inputs, i.e., by utilizing noise-induced synchronization. Various synchronization 
patterns, including fully synchronized and clustered states, can be realized by using linear filters 
that generate appropriate common noisy signals from given noise. The optimal linear filter can 
be determined from the linear phase response property of the oscillators and the power spectrum 
of the given noise. The validity of the proposed method is confirmed by numerical simulations. 


Introduction. — Various nonlinear dynamical sys¬ 
tems tend to synchronize when driven by a common noisy 
input. This phenomenon, called noise-induced synchro¬ 
nization, is observed in many systems, for example, in 
neurons [IHU, electric circuits 0, electronic devices [6], 
microbial cells [7], lasers [5], and chaotic dynamical sys¬ 
tems mm- It has been clarified that noise-induced syn¬ 
chronization has quite a different mechanism from phase 
locking to periodic forcing, i.e., the oscillators are not en¬ 
trained by the input but still exhibit mutual synchroniza¬ 
tion, characterized by coherent distributions of the phase 
differences. Analytical investigations of this phenomenon 
for limit-cycle oscillators can be performed by using the 
phase reduction method m , and its properties have been 
widely studied in the last decade 

In contrast to phase locking that requires periodic forc¬ 
ing whose frequency is close to rational multiples of the 
natural frequency of the oscillators, noise-induced syn¬ 
chronization can occur even for white noise [T2HT81I20U24I . 
Thus, noise-induced synchronization may more easily be 
realized than phase locking, because common environmen¬ 
tal noise is ubiquitous in nature. Indeed, it is conjectured 
that some plants utilize common environmental noise for 
synchronization to realize biological functions [25j[26]- In 
such biological systems, some kind of filtering mechanisms 
for the environmental noise may also exist to improve 
noise-induced synchronization. Such filtering mechanisms, 
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if any, would also be useful in practical applications such 
as noise-induced synchronization of sensor networks |29| . 
The aim of this letter is to provide a theoretical basis for 
optimizing noise-induced synchronization by filtering the 
given noise. 

In our previous works [161123]! we developed quantita¬ 
tive theories that predict global statistical properties of 
the noise-induced synchronization, such as formation of 
various synchronization patterns and fluctuations around 
these patterns, not only whether the oscillators synchro¬ 
nize with the others. The synchronization patterns are 
characterized by the probability density function (PDF) 
of the phase differences between the oscillators, which can 
be calculated from the phase response property of the os¬ 
cillators and the statistical properties of the noisy inputs. 

In this letter, on the basis of ref. [21], we propose 
a method for designing and controlling various noise- 
induced synchronization patterns of limit-cycle oscilla¬ 
tors, including the fully synchronized and clustered states. 
Since the synchronization pattern of the oscillators de¬ 
pends on the statistical properties of the noisy input, we 
can design the synchronization pattern by optimizing the 
noisy input so that some objective function, e.g., the de¬ 
gree of synchronization, is maximized. We develop an op¬ 
timization method for the noisy input, namely, for a linear 
filter that transforms given noise into an appropriate noisy 
input so that the desired synchronization pattern is real¬ 
ized. The validity of the proposed method is confirmed by 


P-1 







W. Kurebayashi et al. 


numerical simulations. 

Model. — We consider an ensemble of N uncou¬ 
pled identical limit-cycle oscillators subjected to corre¬ 
lated noise and independent noise, described by the fol¬ 
lowing Langevin equations: 

X 3 {t) = F(Xj) + eG(Xj)[Ij(t) + <,•(*)], (1) 

for j = 1 Here, Xj(t) £ R 71 is the state of the 

oscillator j at time t, F(Xj) £ R 71 is a vector field rep¬ 
resenting the oscillator dynamics, G(Xj) £ R" repre¬ 
sents the coupling of the oscillator to the noisy inputs, 
Ij ( t ) £ R is the correlated noise, (j (i) £ R is the indepen¬ 
dent noise that drives each oscillator independently, and 
e (< 1) is a small parameter that controls the strength 
of the noisy inputs. We assume that eq. © possesses 
a stable limit-cycle orbit Xo(f) with period T and fre¬ 
quency ui := 2i t/T when e = 0. The correlated noise is 
generated from two given noisy signals by a linear filter 
as Ij(t) = f * (£(t) + r]j(t)), where the noise £(f) £ R is 
common to all oscillators, rjj(t) £ R is independently ap¬ 
plied to each oscillator, /(r) £ 1 is a filter function that 
transforms the given noise to appropriate noise for real¬ 
izing desired synchronization patterns, and the star (*) 
denotes convolution / * a(t) = f+°° f(r)a(t — t)<1t. We 
introduced two independent noise terms rjj (t) and Q (t) to 
take into account the effect of external disturbances before 
and after filtering. 

For example, the oscillators described by eq. © can be 
regarded as spiking neurons receiving artificial injection 
currents as in ref. j3j. In this case, the filtered noise Ij(t) 
represents the injected current to each neuron, and we may 
suppose r]j(t.) = 0. The filter f(r) is used for generating 
an appropriate injection current, and Q(t) is independent 
noise inherent in each neuron, e.g., channel or synaptic 
noise. In the situation of ref. [29] where noise-induced 
synchronization of wireless sensor networks is considered, 
each oscillator described by eq. © corresponds to each 
sensor node. The sensor node measures a noisy environ¬ 
mental signal £(t) + rjj(t), and the filter f(r) implemented 
on each sensor node transforms the signal into an appro¬ 
priate noisy input Ij[t) that induces synchronization of 
the sensor nodes. We may suppose Q ( t ) = 0 in this case. 

We assume that £(t), and (j(t) are mutually inde¬ 

pendent zero-mean Gaussian noise, i.e., (£(t)) = ( r)j(t )} = 
(&(*)> = 0 and - r)) = (i{t)Q(t - r)} = 

(')) = = Ofor 

any j, k , and t (j ^ fc), where (•) denotes the ensemble av¬ 
erage over realizations of £ (t ), r/j (t ), and For simplic¬ 

ity, we assume that the statistical properties of r]j(t) and 
(j(t) do not depend on the oscillator index j. Their power 
spectra are given by P^( fi) := e -7nr (£(t)£(t — r))dr, 

Pr,(n) := J+™ - r))d,T, and P C (D) := 

e~ znT (Cj(t)Q(t—T))dr. We also define the amplitude 
response of the filter /(r) as A(fl) := | e _inT /(r)dr|. 


By the phase reduction method m, we can reduce the 
high-dimensional oscillator dynamics described by eq. © 
to a one-dimensional phase equation for small e, 

dj=uj + eZ(9j)[f * (&(t) + r)j{t )) + Q{t)} 

+ e 2 v(0j) + 0(e 3 ), ( 2 ) 

where 9j(t) £ [0, 27t) is the phase of the oscillator j, 
Z(9j) is a sensitivity function that characterizes the re¬ 
sponse of the oscillator phase to noisy inputs, and v(9j) 
represents the effect of amplitude relaxation dynamics of 
stochastic limit-cycle oscillators [2011323 (this term even¬ 
tually vanishes and does not play a role in the follow¬ 
ing argument). The sensitivity function Z{9) is given as 
Z{9) = G T (X)X x 9(X) \ x=Xo(g/u) , where 0(X) is the 
isochron of the limit cycle and V x 9(X) |x=x o (0/^) repre¬ 
sents its gradient on the limit-cycle orbit at phase * El- 

Characterization of synchronization patterns. — 

As discussed in refs. USE], the phase difference between 
two oscillators <pj t k '■= 0j — Ok characterizes the noise- 
induced synchronized state. Since statistical properties 
of r]j(t) and Cj(t) do not depend on the oscillator index j , 
the PDF of the phase difference 4>j t k does not depend on 
the indices j and k. Thus, in the following, we denote the 
phase difference by cf> without the oscillator indices. 

In our previous work ED, we obtained the stationary 
PDF U((j >) of the phase difference (f> by employing effec¬ 
tive white-noise approximation of the phase equations © 
subjected to correlated colored noise and by deriving an 
averaged Fokker-Planck equation for </> from the multivari¬ 
ate Fokker-Planck equation for the phase variables {6j}. 
It turns out that the correlation functions of the noise play 
an important role, and U(cf>) is explicitly given by 

U((j))=*-—- - 1 ( 3 ) 

u 9(0) - #)+M0) 

where u £ R is a normalization constant determined by 
J’+J 1 ' U{(j))d(j) = 1, and g(<j)) £ R and h{<j>) £ R are corre¬ 
lation functions of the noise terms in eq. 0, defined as 
g{4 ’) = f+~(zm)f * £(t)Z(0{t -T) + 4>)f * £{t - r))dr 
and h(<j)) = J+~(Zm)[f * Vj(t) + Cj(t)} z (0(t - r) + 
4>)[f * + (j(t — r)])dr. In Fourier representation, 

these functions can be written as g{cj>) = )C]£-oo 
and = Y^t^oo hge 1 ^, where the Fourier coefficients 
ge,hg £ I are given by 

9i = M 2 I-4(M| 2 p ? (M, 

hi = | zi\ 2 \A{tu)\ 2 P„{eu) + \zi\ 2 P c (icj), ( 4 ) 

and ze := ^ e~ lM Z(9)d9 is the Fourier coefficient of 
Z(6). From eq. 0, we see that the PDF U(<f>) is sym¬ 
metric about <f> = 0 and has a maximum at (j) = 0. For 
example, a PDF U(<j>) with a single peak at <f> = 0 repre¬ 
sents the synchronized state of the oscillators, and U((j>) 
with k peaks represents the fc-clustered state. Thus, when 
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the amplitude response A(f2) of the filter /(r) enhances 
the /c-th mode of the correlation function g{<j>) in eq. ©, 
fc-clustered distribution is emphasized in the PDF [/(</>). 

Design of synchronization patterns. — Equa¬ 
tion 0 indicates that we can design the stationary PDF 
of the phase difference i.e., the synchronization pat¬ 

tern, by varying the correlation function g((f>). Therefore, 
given the power spectra P v (fl), Pq(£l), and the 

sensitivity function Z(0), we can try to find an optimal 
filter /(t) that gives g(<f>). In this study, rather than ex¬ 
plicitly specifying the precise PDF of the oscillators as 
the target, we aim to maximize its statistical property, 
e.g., the degree of synchronization or clustering. This is 
because such macroscopic properties, rather than precise 
functional forms of U((/>), are relevant in practical applica¬ 
tions. Note also that we cannot generate arbitrary PDFs 
but only optimize the PDF given in the form of eq. ©• 

To design the optimal filter /(r), we introduce the fol¬ 
lowing objective functional R{f} of /(r) characterizing 
the statistical property of U{(j)) as a measure for choosing 
an optimal synchronization pattern: 

R{f}= I*” U(<f)-, f)q{<j>)d(j), ( 5 ) 

J —7T 

where we explicitly show the dependence of /) on /. 
The function q((j>) determines what statistical property we 
focus on. We try to design synchronization patterns with 
desired statistical properties by choosing appropriate q(<j>). 

In the numerical simulations given below, we will use 
the following functions for designing the synchronization 
patterns: = cos <j>, q 2 { 4 >) = <5(</>), <Z 3 (</>) = cos3 <j>, 

and qi{q J) = cos2 (j), where $(•) denotes the Dirac delta 
function. When we use qi(<fi), the objective functional 
R{f} corresponds to the order parameter introduced in 
ref. (18], which characterizes the degree of noise-induced 
synchronization. When we use q2 the objective func¬ 
tional R{f} corresponds to the maximum of the PDF U (4>) 
at <f> = 0, which also characterizes the degree of synchro¬ 
nization, but in a more strict way, i.e., it counts only the 
oscillator pairs with exactly zero phase difference. When 
we use <?3 (</>), the objective functional R{f} characterizes 
three-clustered states, in which three synchronized sub¬ 
groups of oscillators are formed. Similarly, q±{<j)) charac¬ 
terizes two-clustered states. 

In the following, we represent q((j>) as a Fourier series, 
q{<j>) = qo + 2 it cos l(j), where the coefficient g<el 
represents the weight of the £-th Fourier mode. Expanding 
U{(j)) as U{<j>) = ut cosicj), eq. © can be written 

as R{f} = YULi itut- By finding optimal ui for given qe, 
we can obtain a PDF U((/>) and a filter /(r) that maximizes 
the objective functional R{f}. 

Optimization of the filter. — By maximizing the 
objective functional R{f}, we seek for the optimal filter 
/(t). However, unconstrained maximization of R{f} of¬ 
ten leads to divergent /(r). We also need to take into ac¬ 
count that our present theory is not valid for strong noisy 


inputs, because the phase reduction method requires the 
input given to the oscillators to be sufficiently weak m- 
Thus, we should introduce some constraint on the filter 
function /(t). 

In this study, we formulate the constrained optimization 
problem of the objective functional R{f} as follows: 

maximize R{f}, (6) 

subject to a 2 := {Ijit) 2 ) = C, (7) 

where the condition eq. © constrains the variance a 2 of 
the filtered noise Ij(t) to be a constant C. Using the power 
spectra P^(Q) and P V (Q) and the amplitude response A(f2) 
of the filter /(r), the variance a 2 can be written as 

/ +oo 

|H(D)| 2 [P c (H) + P^(H)]dD. (8) 

-OO 

By solving the optimization problem described by eqs. © 
and 0, we can, in principle, obtain the optimal filter /(r) 
for maximizing the objective functional R{f}. 

Actually, we should also take into account that the op¬ 
timal solution of eqs. © and 0 may not be implemented 
in practice. The optimal amplitude response A(f2) ob¬ 
tained as above often has delta peaks at f 1 = iu (£ G Z), 
i.e., |A(D)| 2 = aoS(fl) + Y^T= l a t[8(& ~ + ^ w )] 

(ae is some coefficient and w is the natural frequency 
of the oscillator), because the PDF [/(</>) depends only 
on the harmonic components P^(0), Ppw), P{(2w),... and 
P v (0), P v (lo), P v (2lo), ... of the noise (see eqs. © and ©). 
Such a delta-peaked amplitude response A(f2) corresponds 
to a physically unrealistic filter that extracts only purely 
harmonic components from the noise, which leads to phase 
locking rather than noise-induced synchronization of the 
oscillators. Besides, such singular A(SY) cannot be realized 
in practical implementation of the linear filter /(r). 

To overcome this problem, we restrict the class of A(fl) 
and further assume that the square of the amplitude re¬ 
sponse is expressed as a finite sum of narrow-band basis 
functions as 

m 

|I(D)| 2 := ^ c w W(n-£uj), (9) 

£=—m 

where ci {£ = 0,1,2,...) is a weight coefficient, W(f2) 
represents a narrow-band basis function prespecified be¬ 
fore the optimization process, e.g., a Gaussian function, 
and m is the maximum wavenumber of the filter. We as¬ 
sume that the basis function W (f2) is localized in the range 
|$7| < u>, i.e., W($7) « 0 holds for |f2| > w. The parameter 
m should be sufficiently large to obtain a good filter. The 
restricted amplitude response A(f2) is experimentally fea¬ 
sible, because W(D — £oj) in eq. © can be implemented by 
a band-pass filter that passes frequencies around D = £oj. 

We introduce a new parameter /3 = [/3 0 , /3i, ■ • ■, /3 m ] T 6 
R m+1 as p 0 = \Z\b 0 co\ and /3 e = ^2\bice\ for t = 1,..., m, 
where bi (£ = 0,1, 2,...) is defined as b e = — 
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£o;)[P{(fi)+P,,(f2)]df2. The optimization problem © and 
GD with the above restriction can then be expressed as 

/ + 7T 

U(4>; (3)q(4>)d(t>, (10) 

-7T 

m— 1 

subject to er 2 ss ||/3 || 2 = ^ \/3 e \ 2 = C, ( 11 ) 

1=0 

where U(4>\ (3) is the PDF of the phase difference <fi ob¬ 
tained by plugging eq. ® into eq. d4]), and eq. m follows 
from eqs. © and © and the definition of (3. Thus, we 
can employ (3 as a design parameter of the optimization 
problem described by eqs. m and CD- 

Optimization algorithm. — To solve the optimiza¬ 
tion problem given by eqs. m and m we use the gra¬ 
dient descent algorithm. We randomly choose an initial 
value /3 ( °) and iteratively calculate (3^ for j > 1 , where 
(3^ is the design parameter f3 at j-th iteration. At each 
iteration, we update (3^ as 

pU+1) = pU) + Q v pR{(3 u) ) (12) 

and normalize it as /3^ +1 ) = / 1 |U so that 

/ 3 td satisfies the constraint (fill) . Here, a is a constant that 
controls the step size, and X7pR([3) £ R m+1 represents the 
gradient of R((3) with respect to /3, i.e., V pR{(3) = [ d gg^ i 
MMlT 

• " ’ dpm J ■ 

For simplicity of notation, we define u((j)) := l/[g(0) — 
g((f>) + h(0)}. The normalization constant is given by u = 

u(4>)d(t>. Then, the gradient ^ in eq. (fl2l) can 

be expressed as - = ftl u {<t>)§fc] q (4>) d( l>, 

and the gradients (du{(j))/df3() and ( du/d/3 ^) can be calcu¬ 
lated from eqs. ©, ® and m and the definition of (3 as 

^ 7 ■ = 2« 2 (</>) TJk=o W{ku-£uj) [P € (few) (cos k(j> - 

i)--p,(M], 

Because W(£i) is localized in the range |H| < 

ui, the above expression can be simplified as follows: 

= 2n 2 ^)W(0)^f [P£(M(cos^- 1 ) -P ? (M], 

which reduces the computational cost of the optimization 
process. Using the optimized (3 , we can obtain the opti¬ 
mal amplitude response A(fl) as the square root of eq. ©, 

1/3 I 2 

whose coefficients Co,... ,c m are given by Co = °l = 

for £ = 1,2, Thus, using, e.g., the least- 

squares method E5, we can calculate the optimal filter 
/(t) from the amplitude response A(fl). 

Numerical simulations. — To confirm the valid¬ 
ity of our method, we performed numerical simulations 
using several examples of the objective functionals. In 
the first example, we use the FitzHugh-Nagumo (FHN) 
model of a periodically firing neuron. This model has a 
two-dimensional state variable {u,v), which obeys v(t) = 
v — v 3 /3 — u + Io + I(t) and u(t) = p(v + c — du ), /z = 0.08, 


c = 0.7, d = 0.8, and Jo = 0.875. The frequency of the 
oscillation is approximately w = 0.173 and the noisy input 
Ij(t) is given to v(t). The sensitivity function Z(9) to I{t) 
and its Fourier coefficients are shown in figs. E (a) and (b). 

As the noisy inputs, we use the Ornstein-Uhlenbeck 
noise, whose power spectra are given by P^{^l) = 

P OU (P0.5), P„( Q) = P c (fi) = 0.1 • Pou(n;0.5), wit h 
P ou (^; 7 ) := 7 2 /( 7 2 + ^ 2 ) [fig- E (c)]. As the basis func¬ 
tion IF(fl), we employ a rectangular function, IF(fl) = 
1 (|fl| < w/2), 0 (otherwise). The other parameters are set 
as follows: the coupling strength to the noise is e = 0 . 01 , 
the variance of the filtered noise is C = 10, and the max¬ 
imum wavenumber of the filter is m = 5 (because | ^ | is 
small when i > 6 ). The parameter a used for the gradient 
descent is a = 0.5. Note that the gradient descent algo¬ 
rithm finds only a local optimum and does not guarantee 
global optimality. In order to obtain the global optimum, 
the algorithm should be repeated from sufficiently many 
initial states (3^. 

Figure [2] shows the numerical results for the FHN model. 
For the functions qi{<j>) and ( 72 (</>) defined previously, syn¬ 
chronized states are successfully formed [figs. [ 2 ] (b) and 
(d)]. For the sinusoidal qi(4>), we obtain a filter that em¬ 
phasizes only the first Fourier mode [fig. [2] (a)] (see eq. ©), 
which results in a bell-shaped PDF U{(j>) [fig. [2] (b)]. For 
the delta-shaped q 2 (^), in contrast, we obtain a nontrivial 
filter that consists of multiple modes [fig. [2] (c)]. In this 
case, the PDF U(<p) has a sharper peak than that for qi(4>) 
[fig. [2] (d)] and a more precisely synchronized state is real¬ 
ized. Note that the high-frequency components of A(f2) in 
fig. [ 2 ] (c), which are stronger than the low-frequency com¬ 
ponents, do not significantly affect the statistical property 
of Ij(t), because P(£w) is sufficiently small for large £. 
Therefore, we can safely neglect the high-frequency com¬ 
ponents of A(Q) whose wavenumbers are larger than m. 
For the function <73 (</>), we obtain a filter that emphasizes 
only the third Fourier mode [fig. [2] (e)], which yields a 
three-clustered state as expected [fig. |2](f)]. Note that we 
cannot form a two-clustered state in the FHN model, be¬ 
cause the phase response property of this model has odd 
symmetry, i.e., the second Fourier coefficient 1 2 : 2 1 of Z(6) 
is vanishingly small as shown in fig. E (b). 

In the second experiment, we use the Hodgkin-Huxley 
(HH) model [34] . It also models periodic firing of a neu¬ 
ron, but it has more realistic, higher-dimensional dynam¬ 
ics without odd symmetry, in contrast to the FHN model. 
We apply the noisy input Ij(t) as well as a constant in¬ 
put Iq = 10 to the V variable (i.e., membrane potential) 
of the HH model. The oscillation frequency is approx¬ 
imately ui = 0.438, and the sensitivity function and its 
Fourier coefficients are shown in figs. [3] (a) and (b). In 
addition to qi(</>) and <72 (</>), we use < 74 ( 1 /)) for optimization 
with the aim of forming two-clustered states. The power 
spectra Ppfi), P V (Q) and PpH), the parameters C and 
a, and the basis function W(VL) are the same as before. 
The noise intensity is e = 0.1 and the maximum wavenum- 
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Fig. 1: FitzHugh-Nagumo model, (a) Sensitivity function 
Z(8), (b) Fourier coefficients \ze\, and (c) power spectra of the 
noisy inputs P^(fl), P V (Q) and 
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Fig. 2: FitzHugh-Nagumo model, [(a), (c) and (e)] Amplitude 
response A(fl) of the optimal filter designed by the proposed 
method and [(b), (d) and (f)[ probability density function U(rf>) 
of the phase difference <j> for [(a) and (b)[ [(c) and (d)[ 

52 and [(e) and (f)[ ® (<(>). The insets display snapshots of 
the oscillators in the v-u plane. 
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Fig. 3: Hodgkin-Huxley model, (a) Sensitivity function Z(6) 
and (b) Fourier coefficients. 


ber of the filter is to = 4 (because |z^| almost vanishes at 

£ = 5). 

Figure 0] shows the numerical results for the HH model. 
Synchronized states are successfully formed for qi (</>) and 
92 (0) as shown in figs. 0] (b) and (d). When we use the 
delta-shaped 92 (</>), we obtain a nontrivial filter consist¬ 
ing of multiple modes [figs.0](c)] and the PDF U{<j>) has a 
sharper peak than the case with the sinusoidal q\ (</>) [figs. 0] 
(d)]. In contrast to the FHN model, we can realize a two- 
clustered state as shown in figs. 0] (f), because Z(8) of the 
HH model has a sufficiently large second Fourier compo¬ 
nent as shown in fig. 0] (b). Note that the realizability of 
a particular state is determined by the sensitivity func¬ 
tion Z(9) that characterizes the phase response property 
of the driven oscillator, rather than by the dimensionality 
or complexity of the oscillator model. 

Summary and discussion. — We have proposed 
a method for designing and controlling various noise- 
induced synchronization patterns by filtering the input 
noise, including the synchronized and clustered states. By 
numerical simulations, the validity of the method has been 
confirmed for two types of limit-cycle oscillators. These re¬ 
sults will provide a theoretical basis for optimizing noise- 
induced synchronization by filtering the input noise. 

Though some previous works [27828] proposed optimiza¬ 
tion methods for the phase response property of the os¬ 
cillator to enhance noise-induced synchronization, those 
works considered only the Lyapunov exponent of the phase 
(i.e., the exponential decay rate of the small phase differ¬ 
ence between two oscillators), so that they could not fully 
characterize the synchronized states and could not be used 
to design various synchronization patterns as described 
in this letter. More importantly, in contrast to previous 
works [271128] that gave the optimal phase response prop¬ 
erty of the oscillator, our present study provides a method 
to generate optimal noisy inputs to the oscillator, which 
can be implemented much more easily than designing the 
oscillator response. Thus, our method can be useful in 
various real-world applications, e.g., energy-efficient syn¬ 
chronization control in wireless sensor networks [29] . 

Finally, though we have not considered the effect of dif¬ 
ferences in the natural frequency of the oscillators DZI22I 
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Fig. 4: Hodgkin-Huxley model, [(a), (c) and (e)] Amplitude 
response A(fi) of the optimal filter designed by the proposed 
method and [(b), (d) and (f)[ probability density function U(cp) 
of the phase difference (f) f° r [( a ) an d (b)[ qi {(j>), [(c) and (d)[ 
52 ( 0 ) and [(e) and (f)[ q 4 .((f>). The insets display snapshots of 
the oscillators in the V-m plane, where m is a channel vari¬ 
able [34] . 


[Ml in this letter, it is often significant in practical appli¬ 
cations. Extension of the present method to non-identical 
oscillators will be an important future work. 
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